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Abstract 

In this paper we prove that the partition function in the random matrix model with 
external source is a KP r function. 



1 Introduction 

Let A G H"^" be an n X n Hermitian matrix, and dfi{x) = w{x)dx be a measure on M 
with all moments finite. Then we define the partition function 

Zn{A) = I e^<^^^^Uf,{M), (1) 

JMeH"X" 

where djji{M) denote a unitary invariant measure on H"^" such that if eigenvalues of 
M are Ai, . . . , A„, 

n 

d^{M) = J\^w{Xi)dM. (2) 

i=l 

Due to the unitary-invariance, we assume A = diag(oi, . . . ,a„), Oj G M without loss 
of generality. We consider Zn{A) as a function of eigenvalues of A, and find a KP r 
function property of it. 

Zn{A) arises in the random matrix model with external source [7], [8], [23], |24j . 
Let A G H"^" be an n X n Hermitian matrix, and V{x) be a function defined on M, 
such that e~^(^) decays sufficiently fast. We consider the ensemble of n x n Hermitian 
matrices with the probability density function 

p( Af) = ^(V(M)-AM) 
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where the normalization constant Z^(^;M), also called the partition function, is de- 
fined as 

ZX{A;R)= [ e-^(^W-^*^)fiM. (4) 

The ensemble is called the random matrix model with external source, and we can 
easily identify ZX{A;R) as the Zn{A) in ^ with wix) = e~^'^''\ In the study of gap 
probability [20] , [SI EJ E] , [2] , we need to consider also 

ZX{A;E)=! e-T'-(V^W-^*^)dM, (5) 

all eigenvalues of M G E 

such that the gap probability that all eigenvalues of M are in E, a subset of M, 
is ZX{A;E)/ZX{A;M.). Similarly we identify ZX{A;E) as the Zn{A) with w{x) = 
e~^^^^XE{x), where xe is the indicator function. 

The random matrix model with external source was introduced in the 1990's [TJ [8] 
as a generalization of the standard matrix model, i.e. the A = case, which is first 
proposed by Wigner in the 1950's [16]. However, a special form of the random matrix 
model with external source has been studied by statisticians since the 1920's, under the 
name of Wishart ensemble, one of the most important models in multivariate statistics 

Consider N independent, identically distributed complex samples xi, . . . ,X7v, all of 
which are n x 1 column vectors, and we further assume that the sample vectors Xj's are 
Gaussian with mean and covariance matrix S, which is a fixed nxn positively defined 
Hermitian matrix. If we put Xj's into an n x rectangular matrix X = (xi : • • • : x^r), 
then the sample covariance matrix 5 = j^XX'^ is an n x n positively defined Hermitian 
matrix matrix. If we assume N > n, then the probability density function of 5 is 

PiS) = ^^^e-^'^^^^-''\detSr-, (6) 



is 



where the normalization constant Z^^(S), analogous to the partition function in (jH), 
ZW;^(E) = | e-^^^^-'''\detSf'-dS, (7) 

S is positively defined 

which is the A^S~^) with w{x) = x'^~"'X[o,oo){x). To study the distribution of the 
eigenvalues of S, we also need partition functions like the Z^ {A; E) in See e.g. [1], 

mm- 

The main result in this paper is that the partition function Zn{A) is a KP r function. 
To make the statement precise, we denote power sums of eigenvalues of A 

1 " 

Pfc = ;^5^af, k = l,2,.... (8) 

1=1 

Since Zn{A) is a symmetric function in oi, . . . , a^, it can be regarded as a function of 
{Pi}i£Z+- We should be cautious that pi's are not independent among one another, 
so that Zn{A) cannot be written as a function of piS in a unique way. However, we 
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will eliminate the ambiguity in the following way. First, we define explicitly a function 
ZniPi,P2, ■ ■ ■) in (j26]) , and then identify 

Zn{A)= Znipi,P2,...) (9) 

Thus the rigorous statement is 

Theorem 1. Zn{pi,p2, • • • ) defined by (j26p is a iTP r function in variables pi,p2, ■ ■ ■ ■ 

Since Z„(pi,p2; • • • ) is a KP r function, it satisfies Hirota bilinear equations [9], 
which are equivalent to a series of PDEs. For example, 

which corresponds to the first nontrivial Hirota bilinear equation. 

In [25] Zinn-Justin proved that the HCIZ integral is a Toda r function. See also 
|26j . This result is closely related to theorem [H since the partition function Zn{A) is 
usually evaluated by the HCIZ formula [7J, although we do not use it in this paper, 
and KP r functions are closely related to Toda r functions. Theorem [1] may also be 
proven by the method in |25j . 

Zn{A) was studied as a multi-component KP r function by Adler and van Moerbeke 
et al. In [1] it is proved that the determinant of a moment matrix for several weights, 
after adding deformation parameters, is a multi-component KP r function. By de 
Bruijn's formula (lemma [3]), we find that Zn{A) is a special case of these determinants 
of moment matrices. Especially, when eigenvalues of A have only two distinct values, 
Zn{A) is a 3-component KP r function, and is detailed in [2]. For the nondegenerate 
case, i.e. eigenvalues ai, . . . , o„ of ^4 are distinct, Zn{A) is a (1 + n)-component KP r 
function. 

The intriguing fact is that our KP r function structure is independent of this multi- 
component KP r function structure. In particular, our set of KP flow parameters {pi} 
is none of these sets of KP flow parameters in [T] . 

Remark 1. With the help of Virasoro constraints, it is shown in ^2] that for quadratic 
V, Z^{A;E) in ([5]) satisfies a PDE, which is derived from one Hirota bilinear equation 
of the multi- component KP r function. We can also find Virasoro constraints for 
Zn{pi,p2, • • • ) as the KP r function and get new PDEs satisfied by Z^{A;E). It will 
be done in a forthcoming paper. 

In section [2] we summarize necessary preliminaries. Then in section [3] we give the 
definition of Zn{pi,P2, ■ ■ ■) and prove theorem [TJ 

2 Boson- Fermion correspondence and KP r func- 
tions 



The definition of KP r functions follows that in [14], and all materials on symmetric 
functions are from [211. 
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KP T functions can be defined through representations of the Heisenberg algebra, 
an infinite dimensional Lie algebra. Over any field K with characteristic 0, such as M 
or C, the Heisenberg algebra H denotes the Lie algebra over K, generated by {/ifcjfcez, 
satisfying 

[hk,hi\ = k6k~i. (11) 

We can construct a representation of H over the so called Boson Fock space, which 
is K[pi,p2, ■ ■ ■], the space of polynomials with infinitely many variables, hk's {k < 0) 
act as multiplication operators, and h^s {k > 1) act as derivations: 



hk 



—kp-k for k < —1, 

1 for A; = 0, (12) 

^ for A; > 1. 



This representation of H is called the Boson representation. 

On the other hand, H has another representation over the so called Fermion Fock 
space. To define the Fermion Fock space, we take an infinite dimensional vector space 
V with basis {vijiez- The Fermion Fock space A is composed of semi-infinite forms, 
spanned by the basis 'yio,i_i,i_2,...) which are defined as 

Vio,i_i,i_2,- = Via A Vi_^ A Vi_^ A..., (13) 

with fo,i-i, . . . strictly decreasing and i-k = —k for k sufficiently large. 

We call the semi-infinite form i^o,-i,-2,... = Af_i At;_2 A. . . the vacuum, following 
the physical terminology. Later in this paper, we use form to mean a semi-infinite form, 
unless otherwise claimed. 

Formally, we define the action of hi on V by 

hk{vi) = vi-k, (14) 

and get the induced action on A by 

hk{vio,i_i,i_2,...) = X] ^ ^'-1 ^ '"i-j+i ^ Vi_._k A Vi_._^ A.... (15) 

j>0 

Although the action of hk on V is not consistent with the Lie algebra structure of H, 
the action of hk on A is a representation of H. This is the Fermion representation of 
H. 

We can observe that the Boson representation and the Fermion representation of 
H are equivalent. The correspondence $ between K\pi,p2, . . .] and A is 

$(1) = vacuum, 
^{Pk/k) =/i_jfc (vacuum), (16) 
^(/(Pi,P2/2, . . . )) =/(/i-i, h-2, ■ ■ . )(vacuum), 

where k > 0, and / is a polynomial. Since [hk,hi] = for k,l £ Z~, the polynomial 
of operators /(/i_i, h^2i • • • ) is well defined. Although it is not difficult to check the 
validity of the correspondence the images of monomials on the Boson Fock space 
become messy combinations of the basis of A. It is an interesting question what the 
preimage of t'io,i-i,i_2,.-- i^- The answer is nontrivial, and can be best formulated in 
notations of symmetric functions. 
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Proposition 1. Let iq = Kq, z_i = ki — 1, i_2 = ^^2 ~ 2, ... with kq > ki > K2 > ■ ■ ■ , 
such that {kq, ki, . . .) = k is a partition, then <I>"^(t'io,i_i,i_2,...) = Sk(piiP2; • • • ); where 
Sk is a polynomial in the definition of Schur functions by power sums \j: If we regard 
{pi} as power sums, then the Schur functions s^ satisfies 

Sk = Sk(P1,P2,---)- (17) 

Since every symmetric function can be written uniquely as a polynomial of power 
sums, Sk is well defined. 

To define KP r functions, we need the concept of decomposability of forms. We 
call V £ A decomposable, if and only if 

u = «0 A A u_2 A . . . , (18) 

where Uj's are linear combinations of Vj's. Now we are ready to give the definition: 

Definition 1. f{pi,p2,...) G K[pi,p2, ■ ■ ■] is a KP r function if and only if f = 
^~^{v), where v is a decomposable form. 



3 Proof of theorem [T] 

In order to define Zn{pi,p2, . . . ), we first expand Zn{A) in Schur polynomials. By the 
Weyl integral formula, we have 

J J Ju{n) 

where Cn is a constant, and dU is the Haar measure over U{n). Then we use the 
identity [15j 



e 



Tr{AC/diag(A)C/- 



l{K)<n 

where k's are partitions of k, 1[k) is the length of k, and C^'s are constant multiples 
of Schur polynomials s^, with the normalization 

^C«(X1,...,X„) = {Xi + ---+Xn)^. (21) 



Furthermore, we know [2T], |llj : 



Here the term power sura is defined slightly different from the most common definition in symmetric 
1 

k 



function theory: Pk — \ YlTLi ^ii iii the same style as 
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where H(k) is the hook length product of k, and (n)^ is the Pochhammer symbol: if 
K = {ki, K2, ...,«;/), then 

I Ki 

in)n = llllin-i + j). (23) 
i=ij=i 

Thus if we denote 

G^= j ■■■ j A2(A)s«(Ai, . . . , \n)d^l{\l) . . . dfi{Xn), (24) 

we have 

l(K)<n 

I °° 1 

'-T^Y^Y. 7;;^G«s,(ai,...,a„), 

(25) 



a ^ ^ (n) 

k=0 Khk 
l(K)<n 



and we define 



Zn{Pl,P2,---) = -^Yl Y -A-G^S^{pi,P2,...), (26) 

^"fe=0 Khk ^'^^'^ 

l{K,)<n 



where is the polynomial defined in (jl7p . Then it is clear that Q holds. 

To determine whether Zn{pi,P2, ■ ■ ■) is a KP function, we turn to the Fermion 
representation, and map it as 

1 °° 1 

$(Z„(pi,p2,...)) = 7^V, with V = Y,Y1 T^G^v,, (27) 

fc=o Khk ^""^-^ 

((K)<n 

where for k = (kq, ki, . . . , k;), we denote = Vq+ko, ■■■■ 

To prove that Zn{pi,P2, ■ ■ ■) is a KP r function is equivalent to prove that V is a 
decomposable form. We have a simple criterion (Pliicker relations) for decomposability 
of forms, and first introduce two kinds of linear operators Vi A (v) and i^* (v) on any 
V e A: 

Vi A {vig A Vi_-^ A . . .) =Vi Avif^ A Vi_^ A . . . , (28) 

I i — lVvi,^ A ■ ■ ■ A Vi A Vi . , A . . . if Vi = Vi, , , 
i.x{v^,Av,_,A...)=\\ ^ . (29) 

I U otherwise. 

Here we note that Vi A (v) and ly* (v) are not in A: they are in and A~ respectively. 
A+ is spanned by forms vf- , = fi„ A Vi_, A . . . such that io,i-i, . . . are strictly 
decreasing and i^k = —^ + 1 for k sufficiently large; A~ is spanned by forms v^^ ■ ^ = 
ViQ A Vi_-^ A . . . such that io,i-i, . . . are strictly decreasing and = —k — 1 for /c 
sufficiently large. 

Now we can state the criterion 1141: 
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Proposition 2. v E A is decomposable, if and only if the tensor of forms in A+ (g) A 

oo 

J2 ViA{v)^i^*{v) = 0. (30) 

i=— oo 

For notational simplicity, for a partition k = {kq, . . . , ki), we denote v'^ and v~, 
analogous to v^: 

Vt = Vl+KO,0+Kl,...,-l + l+Kl,-l,..., = V-l+K0,-2+Kl,...,-l-l+Kl,-l-2,...- (31) 

For K a partition and i an integer, we define partitions k + i and k — i. First, 
K = {kq, . . . ,Ki) corresponds to a v ~ G A~, whose subscript {—I + kq, —2 + ki, . . . ,—l — 
1 + —1 — 2,...) is a decreasing sequence of integers. If i G { — I + kqi —2 + ki, — 
1 + Ki,—l — 2, . . .}, then wc say k + i is not well defined; otherwise we can arrange 
elements in {—I+kq, — 2+ki, . . . , —1 — 1+ki, —1 — 2, . . . }U{i} into a decreasing sequence 
which is the subscript of a form v^' G A, and we define n + i = k'. Symmetrically, we 
can define k — i: If i ^ {1 + kq, + ki, . . . , — / + 1 + k;, —I, ■ ■ ■}, then k — i is not well 
defined; otherwise we can arrange {1 + kq, + ni, . . . ,—l + 1 + ki, —I, • • • } \ {i} into a 
decreasing sequence, which is the subscript of a v^" G A, and we define k — i = k". 

Now we consider the tensor of forms 5^^_oo ^ (V) <8) ^v^{V) G A+ A~, and 
compute its coefficients of every t;+ (g) v'^ term. By book-keeping, we get the result 



Lemma 1. For any two partitions a and f3, the coefficient of the f + ® term ofV is 

n(r. R\ \^ sgn+(a,i)sgn'(/3,i) 

f^. {n)a-i{n)i3+i 

a — i and /3 + i are both well defined; 
l{a-i)<n, l{l3+i)<n 

where sgn+(a,z) = (-l)-? if v+ = Vao,a-i,... and i = a-j; sgn-(/3,i) = (-1)-'' if = 
Vbo,b-i,... and b-f < i < 

To simplify c(a,/3), we first prove that 

Lemma 2. For any partitions a and j3, and integer i such that a — i and j3 + i are 
both well defined, we have for any n G Z+ 

{n)a-i{n)p+i = n{n + l)a{n-l)i3. (33) 

Proof. Let a = [uq, ai, . . . , a^), /? = (/3o, /^i, • • • , Pr), and (j = 0, 1, . . . ) 

a_,=|"^-^' + ' (34) 
y-j + l for j > q; 

b^,=i^^-'-' ^""''-'^ (35) 
I —J — 1 for j > r. 
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We assume that i = a_^. and b^i > i > then we have 

(k-in+a^j-l \ / 00 n+a_j-l \ 

n n n n ^ ' ^^^^ 

j=0 p=n-j J \j=k+lp=n-j+l J 

n n f n f n w A p^) 

i=0 p=n-j I \p=n-l-l I \j=l+lp=n-j-l I 



On the other hand, we have 



Here we take the convention that J^^, = 1 if < p' . 



(k-in+a-j-l \ /n+a_fe-l \ 
n n H n H I n n (^8) 

j=Op=n— j+l / \p=n~k—l I \j=k+lp=n—j+l I 

(I n+fe_j-l \ / oo \ 
n n H n w A (39) 

jr=Op=n-j-l / \j=l+lp=n-j-l I 



Thus we get 



(40) 



(n + i).(n-i), nU--^--i)nr:;r^' 

Notice that a_fc = i, we can verify that 

(n)„,,(n)^+, = ^ (41) 
(n+l)«(n-l)^ ' ^ ^ 

and prove the lemma. □ 

From this result we observe that if both the a — i and f3 + i are well defined, the 
condition max(/(a — i), l{(3 + < n is equivalent to max(/(a) — 1, + 1) < n. Then 
we have for l{a) < n + 1 and l{j3) < n — 1, 

c{a,P)= I -I A X] sgn'^{a,i)sgn'{P,i)Ga-iGf3+i, (42) 

o — i and /3 + i are 
both well defined; 
max(/{a)-l,«(/3)+l)<n 

and later in this section we assume max(/(a) — + 1) < n. We find that for all 

but finitely many i E Z, either Ga-i or G/j+j is not well defined, and we can write 
as a finite summation 



1 

c(a,/3) = , , / 7T- X]^Sn"^("''^-j)^Sn~(/3,a_j)Ga_a_jG'/3+a_^., (43) 

n(n + l)a[n ^.^^ 

where a-j is given by (p^ . and we assume Gj^j^a^j = if /? + a_j is not well defined. 
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Now we recall the determinant al formula for Schur polynomials |22] that if k 
(ko, ■ ■ ■ ,Ki) with / < n, then 



Sk , . . . , Xn) 



n-l+Ko 
n~2+Ki 


n-l+Ko 

• • • 

n-2+Ki 
^2 


■^n 


-1 + Ko 


n-/-l+Ki 


n-/-l+Ki 
X 2 ... 




-l + K; 
-i-2 


1 


1 




1 


A(xi, . . . ,a;„) 



(44) 



Therefore 



•vn-2+Ki \n-2+Ki 
Al . . . A„ 



A(A)2s,(Ai,...,A„) = A(A) 

1 ... 1 

To simplify the integrals in ()42p . we need another formula [TOl 

Lemma 3 (de Bruijn's). For any /o, . . . , /n-i, 5o, • • • , 9n-i G 



(45) 



dXi . . . dXn 

/n-l(Al) /n-l(A2) ••■ 

/n-2(Al) /n-2(A2) •• 

/o(Ai) /o(A2) ... 



Now we denote 



/n-l(An) 
/ft-2(An) 



fi'n-llAl) 5n-l(A2) 
5'n-2(Al) 5'n-2(A2) 



5n-l(An) 
5n-2(An) 



/o(A„) 5o(Ai) 5'o(A2) 

n!det ( / fi{x)gj{x)dx 



0<i,j<n-l 



(46) 



= y xM^(x), (47) 
and by ()45p and lemma El for < j < n we have (a_j is given by ([3^ 



Ga-a., ^ J " J ^(A)^SQ-a_j (Al, . . . , An)d)u(Ai) • . . dAt(A„) 



n! 



Gn-l+a_i Gn-l+a^i+1 



G*}!— 1+ao+n— 1 
G*}!— l+a_i+n— 1 

Gn-l+a_,+n~l 



l+a_„+n— 1 



(48) 
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where " means the entry is deleted. For /9 + a_j, we denote = bb^^h^i,... and similarly 
have that if /? + a_j is well defined and /(/3 + a-j) < n, then 



sgr 



= s 


gn {P,a^j) J 


[■■■J A{Xfsp+a. 


_^-(Ai, . . . , Xn)dn{Xi) 






Gn-l+a^j+1 


Cn— l+a_j+n— 1 




Gn-l+bo 


Gn-l+bo+1 


Gn~l+bo+n-l 


n! 




Gn-l+b-i+1 


Gn~l+b-i+n~l 








G„-l+b-„+2+n-l 



(49) 



Here we notice that the nx n matrix in (I48p is constructed from an (n + 1) x n matrix 
with the (j + l)-th row eliminated, and the first row in the matrix in (|49p is the same as 
the deleted row in the construction of the matrix in (I48p . Thus we have a determinantal 
formula 



^sgn+(a,a_j)sgn {(3,a-.j)Ga~a.^Gp+a.j = 

j=0 

Bq Bi 

Gn-l+bo Gn-l+bo+1 
Gn-l+b-i 



n! 



G'n-l+b_i+l 



G 



n—l+b-n + 2 Gn—l+b-n + 2 + 1 



Bn-1 
Gn-l+bo+n-1 
Gn-l+b-i+n-1 



G. 



n-l+b-n+2+n—l 



(50) 



where all rows except for the first one are the same as those in the matrix in (|49p . and 
the first row {Bq, Bi, . . . , B^-i) is 



(-Bo, Bi, . . . , Bn-l) = 

n 

^ Sgn+(Q, a-j)Ga-a^j (G ), (51) 

j=0 

Notice that sgn"*"(a, a_j) = (—1)-', we find that (j5ip is equivalent to that Bf^ is the 
determinant of a (n + 1) x (n + 1) matrix (/c = 0, 1, . . . , n — 1): 



Bk 



Gn-l+ao+k Gn~l+ao Gn-l+ao+1 
Gn-l+a-i+k Gn-l+a^i Gn-l+a^i+1 



Gn—l+ao+n~l 
Gn-l+a^i+n-1 



Gn~l+a- 



Gn- 



Gn- 



l+a_„+l 



Gn-l+a_ 



(52) 



Since in the matrix, the first column is identical to the {k + 2)-th column, we get 
Bk = 0. Therefore ^"^q sgn+(a, a_j) sgn~(/?, a„j)Ga_a_^ G'/3+a_^. = 0, which means 
that c(a,/3) = by (03]). 

Summing up these results, we get the conclusion that V is a decomposable form, 
and by (j27p prove theorem [TJ 
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